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1. Introduction 

It is by now well established that react ion- diffusion models can have a rich dynamics 
governed in low dimensions by the fluctuations. Several physical quantities behave 
as power law and the associated critical exponents have some universal properties. 
The renormalization group (RG) method developed in the framework of equilibrium 
statistical physics fl], provides a suitable tool to study such dynamics. To be able to 
approach the problem of reaction-diffusion in term of RG, one has first to go from the 
initial microscopic master equation to a coarsed-grained description. The standard way 
to do this consists in using a Fock space formalism (see the works of Doi ||, Grassberger 
and Scheunert and Peliti Q). One ends up with a model which dynamics is defined by 
the action of a continuous field theory. Among several applications, this method has been 
used to study the two species annihilation reaction problem A+B — > 0; for homogeneous 
initial state with equal densities 0, or unequal densities 0, rigorous predictions of 
Bramson and Lebowitz have been reproduced and novel results obtained. For the 
case, in which the two species are initially spatially separated, Howard and Cardy ||, 
have confirmed scaling arguments developed by Cornell and Droz ||. 

Another interesting family of reaction-diffusion processes is formed by the one- 
species diffusion-annihilation and the diffusion-coagulation models. In 1986, Peliti [ |TC[ 
showed that the coagulation-diffusion model A+A — > A belongs to the same universality 
class than the annihilation- diffusion one A + A — > 0. He also showed that the associated 
field theory is super-renormalizable. Thus, only the coupling constant needs to be 
renormalized. Moreover, this renormalization can be done to all orders in e = d u — d, 
where d u , the upper critical dimension is 2 in this case. Peliti showed that the 
concentration c(t) of the reactant in the long time regime behaves as: 

c(t) oc t~ a , a = ^, (t — > oo). 

However, he made no predictions concerning the amplitude, neglecting the initial 
conditions in its approach. It turns out that the initial conditions may play a very 
important role. This aspect has been taken into account within this formalism only in 



1991 by Ohtsuki [[nj and in 1992 by Friedman et al pf. In 1994, Lee |J gave the first 
complete RG analysis of the annihilation-diffusion model. It turns out that the initial 
conditions show up as a local source into the action. Lee was able to treat this term to all 
orders in perturbation theory. It was shown later (see for the two-species annihilation 
and JT4J for the one-species case) that such infinite resummation is equivalent to a shift 



of the fields in the action by their classical values. 



Meanwhile, in 1993, Droz and Sasvari [15] addressed the problem of both 



annihilation-diffusion and coagulation-diffusion in presence of a source J of particles: 
A + A^®, A + A 1 * A, 0^>A 
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Performing a renormalization procedure, they found that the density of particles obeys 
the scaling law 

7 \ d/(d+2) 

J \ — " \2/(d+2) 



c(t,J,D) = ^-) F[(J/D)W)Dt), (1) 

for sufficiently small value of J. However, the scaling function JF was not computed. 
Making ad hoc assumptions on the asymptotic behavior of T in the limits t — > oo, they 
showed that the stationary particle density 

c-J 1 /*, 6=1 + 1 

a 

was approached with a characteristic relaxation time r given by 

2 



J- A A 



d + 2 



Moreover, considerations on the J — > limit, allow them to reproduce the scaling laws 
postulated phenomenologically by Racz [HJ, namely: 

a5A = l, A+- = l. 

o 

The goal of this paper is two folds. First to provide a complete renormalization 
group analysis of such models by computing not only the exponents but also the scaling 
function T defined by ([I]) in the framework of an e = 2 — d expansion. This is 
done in section 2. Second, to give an exact analytical solution of the one-dimensional 
coagulation-diffusion model with infinite reaction rate and source, by extending to time 
dependent regime an approach developed by Doering and ben-Avraham 17]. These 
exact results are compared with the RG predictions in the limit e = 1 in section 3. 
Remarks and conclusions are given in section 4. 



2. Field theoretical approach and renormalization group analysis 

2.1. The model and the associated field theory 



We shall not derive here in details how one obtains the field theoretical model. 
The interested reader is referred to the original papers of Doi [0, Grassberger and 
Scheunert [[| and Peliti 0, |Kj or to the short reviews presented in |13| and in ||18|| . 
As the number of particles is not conserved, the basic idea is to introduc a Fock space 
representation. The time evolution operator of the problem can then be cast in a path 
integral form which, in the continuous (coarsed grained) limit is characterized by the 
action: 



5* 7 [a, a, J] 



l d x f 
Jo 



dt 



a + •y\aa z + Xa 2 a 2 — J a 



(2) 
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A whole class of models indexed by the parameter 7 G [1,2] is thus defined. For 7=1 
one has the pure coagulation-diffusion model and for 7 = 2 the pure annihilation one. 
For 1 < 7 < 2, both reactions are possible with a given probability depending on 7 
(see for more details). The coupling constant A is related to the reaction rates g 
and k via A = (7 + l)k and 7 = (7 + 2) (7 + with 7 = k/g. The particles diffuse 
in an infinite <i-dimensional space with a diffusion constant D. The above action could 
models two different types of colliding particles in some appropriate limits. First point- 
like particles living on a d-dimensional hypercubic lattice JIB] , second extended particles 
living in a rf-dimensional continuous space JT8[ . 



The time and position dependent fields a and a obey bosonic like commutation 
relations. The field a(x, t) is related to the local particle density, while the auxiliary 
field a(x,t) has no particular physical meaning. 

Within this formalism, correlation functions are expressed by functional integrals. 

G^({x u = JvaVa a(x l ,t 1 ) . . . a(x N , t N ) 

xa(x N+ i,t N+1 ) . . . a(x N+] y, t N+] y) exp(-S a [a, a, J]). (3) 
In particular, the particle density c 7 (x, t) at point x and time t reads: 

c y (x,t) = JvaVaa(x,t) exp(—Sj[a,a, J]) (4) 
As we have shown in a previous work |18| , the correlation functions of different 7 models 



are closely related. In particular, one has 

S y [a,a, J] = 5i [7a, 7 _1 a, 7 J], 
which implies for the concentration 

c 7 (x, t; J) = 7~ 1 c :L (x, t; 7J). (5) 

Accordingly, it suffices to study one particular model belonging to the class to know 
the behaviour of the other members. From now on, we shall study the pure coagulation- 
diffusion model (7 = 1), with an initial state empty, which action is: 

S[a, a, J] — J d d xj*di a(J^- DV 2 ^ja + Xaa 2 + Xa 2 a 2 - J a , (6) 

For the sake of simplicity, we do not write the index 7 = 1 in the following. When 
J = A = 0, one has a free theory (pure diffusion) and the spatial Fourier transform 
of the free propagator is simply G (p,t) = 6{t) exp(— Dp 2 t), where 9(t) is the usual 
Heaviside function. Simple power counting shows that the upper critical dimension 
of action (||) is d u = 2. For d > d u , the quadrivertex Xa 2 a 2 is irrelevant and can be 
discarded. Below d n , the quadrivertex Xd 2 a 2 is relevant and leads to singularities that 
have to be renormalized. At d = d n , this vertex is marginal, and one expects logarithmic 
corrections to the mean-field behaviour. 
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2.2. Mean-field solution 

We first consider the case d > d u = 2, where the behavior is mean-field like. At the 
mean-field level, the equations of motion for a and a are obtained from the action (|5p, 
by the usual saddle point argument and read: 
5S ( d 

and 

5 a \dt 

Assuming that a and a are homogeneous, it follows that as expected a = is a solution 
and (0) becomes 

= -\ a 2 + J 

at 

with the initial condition a\t=o = 0. Thus, the mean-field or classical solution is: 



^= [^--DV 2 )a + Xa 2 + 2Xaa 2 - J = (7) 
5a \ot 

= -f^- + DV 2 )a + 2Xaa + 2Xa 2 a = 0. (8) 



= J± 1-2 «p(-W 
W V A \ l + exp(-2v / JAt)/ 
where the subscript "cl" stands for classical. 

2.3. Renormalization 

Let us now consider the problem below 2 dimensions. A brute force computation of 
the particle density from (f|) leads to divergences. Thus the coupling constant A needs 
to be renormalized. Neither the fields a and a nor the diffusion constant D require 



a renormalization |10|, [L3|], and as a consequence, the particles input rate J does not 
either. 

We define the temporally extended vertex function X(p,t) to be the sum of the 
diagrams shown in figure [I]. These diagrams can be summed to all orders and the 
Laplace-transformed vertex function reads (see fill): 

^' S) = 1 + 2[X/(8nDy/ 2 ]T(le)(s + \D V 2 )^/ 2 (10) 

We define #r, the dimensionless renormalized coupling (or running coupling constant), 
using the minimal subtraction scheme. That is we define 

gji = Z s K~ e -^- (11) 
yR g 2nD V ; 

where k is a normalization point and Z g = 1 + Y^ili a i9R- The are chosen such as to 

exactly cancel the poles of order 1/e appearing in A(p, s). From (0)> we have 

A 



A(p, s) 



k-*X[1 + tA{e)]/{2TiDe)' 
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+ -<xxx +■■■ 

Figure 1. Diagrammatic representation of the temporally extended vertex function 
A(p, ^2 ~ The propagator is represented by a full line. Here is shown the 

diagrammatic expansion for the trivertex Xaa 2 . Similar diagrams can be drawn for 
the quadrivertex Xa 2 a 2 . 



where A(e) = 0(1). With ([IT]), one finds 



A(p,s) k e g R 



2nD Z g + g R /e + g R A(e) ' 

Accordingly to our prescription, we have to choose a x = — 1/e and a { = for i > 1, 
giving Z g = 1 — g R /e (exact to all orders). The j3 function is defined by 

a, \ _ &9R 2 

P[9r) = k-q^ = ~e#R + fe- 
lt is exactly quadratic in g R and has a fixed point given by /3(g R ) = at g R = e. 
The bare coupling may thus be expressed in terms of the renormalized one: 

= 9R t =g R + ^ + ... (12) 
2nD l-g R /g R 9R g R 1 J 

The perturbation theory can then be written as an expansion in powers of g R . 
2.4- Renormalization group equations 

An arbitrary renormalized correlation function G R N ({xi, t%}i +N ) (where the subscript 
R stands for renormalized) is related to its bare expression (3) through 

9b, D, J, «) = G N >*({x h U}^; A, D, J). 



The independence of the bare functions on the normalization scale can be expressed via 
the condition: 

d . „, 

&9R' 

The formal solution (obtained by the method of characteristic) is 

t,}f g R , D, J, k) = t t }^; g R (p), D, J, pn) (13) 

with 

gM = g R (l + 9 -^^A'\ (14) 

V 9r / 
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Note that in the small p limit, g R (p) —> g R . 

We can implement (|T3| ) with a dimensional analysis. The dimensions of the different 
quantities, expressed in term of momentum k and energy E are 

[t] = E~\ [D] = Ek~ 2 , [J] = En d 

[a)=K d , [5] = 1, [G^({x l ,t l }^)] = n Nd , 

Thus, 

G^ax^+^^^D, J,k) = K Nd G^({Kx t ,K 2 DU}? +R ;g R ,l,K- d - 2 J/D,l). (15) 
The combination of (|H|) and flT5|) leads to 

= (p^G^W,, (pKfDti}" g R (p), 1, (p^)- d " 2 J/D, 1). (16) 

We can then use the following strategy to compute the correlation functions: first an 
expansion in power of A is established; then it is converted into an expansion in power of 
g R through (|l|). The singularities in e are eliminated using the renormalization scheme. 
Now, for a correctly renormalized theory, we can rewrite the g R expansion into an e 
expansion using (16) and (0). Indeed, introducing the p dependence through (16) and 
letting p — > 0, <7r(p) — > g R one obtains G R N as an expansion in power of e. 

Up to now p is an arbitrary parameter, and several choices are possible. For example, 
if we choose p such that 

(p*r d - 2 i = !■ ( 17 ) 

the limit p —>■ becomes equivalent to J — > 0: when the source rate is vanishing small, 
the running coupling approaches its fixed point value. 
Another choice is 

{pKfDtt = 1 

and the limit p — > may be exchanged with t\ — > oo. However, from (16), we see that 
(pK)~ d ~ 2 J /D diverges. Thus, with this choice, one should know the behaviour of G^ N 
for arbitrary large value of J. However, for large values of J, the e-expansion breaks 
down (see Appendix A). Accordingly, we shall choose p acccording to condition flUTD in 
what follows. 



2.5. Density calculation 

The density is first calculated using (^), in the framework of a perturbation expansion 
in power of A. At the tree- level, we find out the mean-field result and we can directly 
use the RG equation (no renormalization is needed). The first correction to the classical 
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behaviour is given by the one-loop contribution. The corresponding diagram may be 
calculated using the action obtained from (0) by shifting the field a by its classical value: 



S[q, fj,J} = J d d xj^ dt v(i~ DV 2 + 2Xa cl ^jri + XfjT] 2 + Xfj 2 (a 2 d + 2a cl r] + 
where r\ = a — a c \ and fj = a, with a c \ given by @. We have 



G m {p, t, f) = 6{t - t') exp[-Dp 2 (t - t')} 



'cosh( v / JAt') 



cosh(y JX t) 



Note that because of the initial condition, G vf j is not invariant under time-translation. 
Obtaining divergent expressions for the one-loop corrections, we shall renormalize them 
using the renormalization scheme developed above. We give below a summary of these 
results. 



2.5.1. Tree level Applying the RG formalism as developed above on the mean-field 
equation, we find for J sufficiently small, 

1 /J\ d /( d + 2 ) 
c R (t; 9R ,D,J,.) = - / =^-j 



f 2eM-2V2^(J/D) 2 ^Dt] \ 

{ 1 + 2 exp[-2 v / 2^ [J/Df/^Dt] ) [ + U[)l ' 1 J 

which is universal (independent of g R )- This result is valid for any time t, because we 
only need to tune J to be in the critical domain. 
For small t we find the expected result 

c R (t) = Jt + 0(t 2 ) 

and for long time (t — *■ oo) 

c R (t) = ■j=(p) {1 - 2exp[-2v / 2^ (J/D) 2 ^ Dt]}[l + 0(e)]. (20) 

The steady state value is thus given by 

1 / J\ d /( d + 2 ) 

c « {co) =V57(d) |1 + oWI (21) 

and it is approached exponentially in time as: 

5cB.it) = c K (t) - c R (oo) oc exp[- % (J/D) 2 ^Dt] (t - oo) (22) 



with % = 2 v / 27re [1 + 0(e)]. 
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t t 2 h 




-p 



Figure 2. One-loop diagram for the density, using action (|18|). The double line stands 
for the free propagator and the dot for the vertex Xa^fj 2 . 



2.5.2. One-loop corrections The diagram corresponding to the one- loop correction is 
given in figure U its analytic expression is 



C W(f) 



2JA 



coslr(vJAt) 



coslr ( V J\ h) smh z ( V JX t 



cosh 2 (VJA t 2 ) 

The integral over the momentum gives the factor [87r.D(t 2 — ti)]^ d ^ 2 - We thus obtain 
2JA t 1+ ^ 2 



c (1) (t) 



/ dx 2 / dxi x 1 1+e / 2 
Jo Jo 



{8nD)- d / 2 cosh 2 (VJ\t) 

cosh 2 [v^JA £2(1 — x{)] sinh 2 [-\/jAa;2(l — Xi)] 



x ■ 



cosh (v JX x 2 ) 

In the framework of an e-expansion, we eventually find (treating generalized 
function pPf , see Appendix B): 



c«(t) 



where 



2JX 



t l+e/2 



{tocD)-*/* CO sh 2 (VJAt) Ve 



2sinh(2VJAt) - 2VJAt 



2VJXt 



+ 



JXt) 



JXt 



O 



0(0 = £ / dx 2 lnx 2 sinh 2 (£x 2 ) 
Jo 

- 1 dxi f 1 , /cosh 2 [£ x 2 (l — Xi)} sinh 2 [£ x 2 (l — Xi)] 
1 ax-? ' 



Jo 



Xi Jo 



cosh (£,x 2 ) 



- sinh 2 (£x 2 ) ). 



Putting together the mean-field result and the one-loop correction, we verify that the 
density is indeed divergence free, and using the RG equation (16) we find, for small 
values of J 

1 



pr(*; 9r, d , J, K ) 



+- 



J\d/(d+2) 



t&nh[V2^e(J/D) 2/{d+2 >Dt] 



cj>[V2^ {J/Dfl^Dt} 1 
cosh 2 [ v / 2^(J/^) 2 /( d+2 )^] + 2 



\n[87i(J/D) 2/(d+2) Dt] 



x 



tanhfv'^re (J/D) 2/(d+2) Dt] v ^ 7 ; 



cosh 2 [V2^ ( J/D) 2 /( d + 2 )Dt] , 



+ 0(e^ 
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J \ d/(d+2) 

T[{J/Dfl^Dt\. 



(23) 



We immediately identify the scaling function T defined in (|I|). Taking the limit t — > oo 
we find 



cr(oo; g R , D, J, k) 



1 



J \ d/(d+2) 



l- 6 - 

4 



Te 



-Kt) 



+ 0(e^ 



(24) 



where 7e is the Euler constant (te — 0.5772). Unfortunately, due to the complicated 
form of 0, we are unable to give a more compact form for the asymptotics of T . Notice 
that the empty initial condition implies that limj_> crW = 0. 

Let us consider two particular cases of interest. The first one is the case e = 1 
(d — 1). Then the steady-state density is 

cr(oo) = -^(^) V3 [1 - l(T E - InSrr) + • • ■] (25) 

and it is asymptotically approached (at the tree-level) as: 

5c R (t) oc exp[-V8^(J/D) 2/3 Dt}. (26) 

We shall compare the accuracy of these expressions with exact results in the next section. 
Note however that, in principle, nothing ensures us that the terms we neglected are small. 
The second case is e = (d = 2). The running coupling is given by 



9r(p) 



9r 



1 - £a In p 

which goes to — (hip) -1 when p — > 0. By replacing this expression in our previous 
formula, we find for the steady-state density 



Cr(OO) 



2jy/ 2 



\n(J/n 4 D) 



-1/2 



ohHJ/^D)]- 1 



and 



5c R (t) oc exp[-v / 8^(J/D) 1/2 \n(J/K 4 D)Dt]. 
as anticipated, logarithmic corrections to the mean-field result are obtained. 



(27) 
(28) 



3. Exact results in one dimension 

A large amount of work has been done to solve exactly one-species diffusion reaction 
models in one dimension (see for example In particular, the diffusion- annihilation 

and the diffusion-coagulation reactions have been considered with an input of particles. 
In the diffusion-annihilation case, Racz |22] obtained the steady-state concentration by 
mapping its model to the kinetic Glauber-Ising model P3] . In 1988, Doering and ben- 
Avraham [23]] calculated exactly the time dependent concentration for a simple diffusion- 
coagulation model, using the interparticle distribution function. Since, their method 
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= J 

t=0 



has been generalized to other diffusion-coagulation processes (see for example [fffil ) and 
in particular the steady-state concentration has been obtained [HI] for the diffusion- 
coagulation with an input source of particle. 

In this section we aim at testing the validity of the RG predictions for d — 1. For 
this purpose we shall extend Doering's and ben-Avraham's results and compute the 
time-dependent concentration. 

We consider an infinite chain (our one- dimensional space) initially empty, and we 
allow particles to appear randomly at rate J (per unit time and per unit length). Thus 
initially 

dc{t) 
dt 

where c(t) is as before the particle concentration. The particles diffuse on the line (with 
a diffusion constant D) and when two particles meet, they instantaneously coagulate 
(A + A — > A). Note that this model is the same as the pure diffusion-coagulation process 
of section 2, but with an infinite reaction rate A. Of course, for such a reaction rate, the 
perturbation expansion in power of A is meaningless and one may argue that the two 
models are not equivalent. However by examining the relation between the renormalized 
coupling g R and A, we see (from (0)) that: 

g R = g* R (l + 2rcg* R ^D/X)- 1 ; 

and thus, when A is infinite, g R = g R . We argue that the infinite reaction rate limit may 
be obtained by taking the fixed point coupling limit (#r — > g R ), that is to say by taking 
the A —>■ oo limit after having performed the path integrals. This will be confirmed in 
the following, at least in one dimension. 

In a one-dimensional space, the particles concentration can be related to the 
probability E(x, t) that an interval of length x > is empty at time t, via 

dE(x,t) 
dx 

As shown by Doering and ben-Avraham [17|, E(x, t) has the advantage to obey a closed 
equation of evolution, namely: 

< 29 > 

with the two conditions 

E(0,t) = l, E(oo,t) = 0. (30) 

From this equation, one can immediately obtain the steady-state, by setting the left 
hand side to zero. One then recognizes the Airy equation, whose solution is (taking into 
account conditions (p0|)) 

A\Uj/2D) l ^x 
E ^°° ) = Ai(0^ 



c(t) 
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where Ai(z) is the Airy's function (see PEfl ). As a consequence, the asymptotic 
concentration reads 

Ai'(o)/ jy/3 

c( °° ) = -AiWfeJ (31) 

Ai'(z) is the first derivative of Ai(z). Note that Ai'(O) < 0. 

Before comparing (31) with the RG results, we shall compute the time dependent 
part of the concentration. For this purpose, we shall solve ( p9[) using the Laplace 
transform E(x,s) defined by E(x,s) = J^dt e~ st E(x,t). Equation (|29| ) becomes 

2D 92E ^ S) - (Jx + s)E{x, s) + E(x, 0) = 0, (32) 

where E(x, 0) is the initial condition (for an empty system, E(x, 0) = 1). Equation ([32]) 
is an inhomogeneous second order ordinary differential equation. Its general solution 
is the sum of the homogeneous solution E^x, s) and a particular solution. The 
homogeneous solution is 



Eh(x, s) 



\2D 



[ai(s)Ai(z) + a 2 (s)Bi(z)\, 



where z = (J/2D)- 2 / 3 (Jx + s)/2D, ai {s) and ct2{s) are two unknown functions of s and 
Bi(z) is the second Airy's function ||26|| . 

Writing E(x,s) = {J/2D)~ 2 ^A(z,s), © becomes 

d 2 A(z,s) , 1 

■ - zA(z.s) = , 

v ; 2D 



dz 2 



dv [Ai(v)Bi(z) - Ai(z)Bi(u)]. 



for which a well known solution is 

A(z, s) = , 

v ' ; 2D Jo 

The two boundary conditions (|30p permit us to determine the two unknown functions 
ai(s) and ^(s). We eventually find for the general solution of 



E(x, s) 



7T / J 

2D\2D 



2/3 1 

Ai(e + a 
Ai(a) 



(Ai(e 
1 



cr 



dv Bi(f ) + Bi(f + a) 
V 



dv Ai(v) 



ira 



+ Bi(a)(J°dvAi(v)--^ 



(33) 



where £ = (J/2D) 1 ^x and a = (J/2D)- 2 ^s/2D. 

The probability E(x,t) is then obtained by Laplace-inverting (33). For t > 0, we 
only have to care for the poles of E(x, s). They are located at a = and at a — a n , 
n = 1, 2, 3, . . ., where a n is the n-th zero of Ai(x) (a n < 0). We finally obtain 



E(x,t) 



Ai(0) 



£ 

71=1 



Aj(g + Or, 

Ai'(a„) 



- 1 








1- 7rBi(a n )| 


v./) 


- 





df Ai(v) 



-|<1ti|t 



with r = 2Dt(J/2D) 2 / 3 and 

4 t) = (±\*( Ai'(O) 



h 7rBi(a n ) 



df Ai(v) 



-|a„ 



\2DJ I Ai(0) n=1 

We are now in position to compare these results with the RG results for e = 1 
the steady-state, the RG gives (up to the one-loop corrections) (EH) 



Cr(OO) 



0.531 ' 



1/3 
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(34) 
For 

(35) 



(36) 



(for small J) while (putting t = oo inside (l34j )) the exact solution gives: 

. Ai'(o)/ jy/ 3 n , 8 /</\ 1/3 

(for arbitrary J). Surprisingly, the difference is only of the order of 10%. 

The comparison for the approach to the steady-state is less convincing, mainly due to 
the fact that we do not know the one-loop corrections. The RG gives, from equation (f26|) 

r i/2\ 1 /2/J\ 1 / 3 

Sc R (t) = - 1-I( 7E -M27r) ' ' ' 



7T 



D 



— exp[-V8^(J/D) 2/3 Dt] 



106l i 



1/3 



exp[-5.01(J/D) 2/3 L>t], 



(37) 



whereas the exact result is 



Sc(t) 



<ii 



dv Ai(v) 



h 7rBi(ai 

ai 

-1.10 f— J exp[-2.95(J/D) 2/3 J Dt] 



f- 

\2D 



1/3 



|oi|t 



(3? 



(with ai ~ —2.33). Both amplitudes in front of the exponential are in good agreement 
(because we used the one- loop result of the steady-state). However, the tree- level 
amplitude into the exponential is quite different (of almost 70%) from the one given 
by the exact theory. The inclusion of the one-loop correction should lead to a better 
agreement. 

Note that the exact results are valid without any restriction on J, on the contrary to 
the RG results, which apply only for small J. This restriction was introduced to ensure 
that (7r to be inside the critical domain (i.e. near ppj. In view of the predictions of 
the exact theory, it turns out that this restriction over J is unnecessary. The coupling 
stays inside the critical domain for any value of J. This justifies why when A — > oo, 
9R -> 9r- 



4. Discussion and conclusive remarks 



We have shown that the RG is a suitable formalism to compute the density in a wide 
class of diffusion-reaction model, with an input of particle. In particular, we easily 
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calculated the critical exponents in arbitrary dimension. However, the computation of 
the universal scaling function is generally much more difficult (see for example (23)) and 
can only be achieved within a power expansion in e = 2 — d. However, the comparison 
of the first order results for e = 1 are not to far from the results obtained exactly in one 
dimension by a different approach. 

One of the main interest of the RG method is that, contrary to most unidimensional 
exact approach, the RG approach is not restricted to the computation of the 
particle density but higher correlation functions can also be calculated (although the 
computation can become quite involved, see [|13|] ). 

Another advantage of RG approach, is the fact that the properties of a whole class of 
model can be simply related. This allows us using (|j|), and knowing the particle density 
for the diffusion-coagulation model, to obtain the density of any mixed annihilation- 
coagulation process. In particular, for the pure diffusion-annihilation model (a = 2) in 
one dimension, one recovers the steady-state density previously calculated by Racz |22|] . 

The present work can be generalized in several directions. Generalization to mA — > 
(m > 2) reactions with a source of particles is also straightforward. Indeed, following 
Lee, it turns out that apart the lowering of the critical dimension to d u = 2/ (m — 1), only 
minor changes occur (for example the amplitudes become m-dependent). This is due 
to the fact that in our renormalization scheme only the fixed point of the coupling 
is modified. The structure of the equations remain unchanged and in particular, the 
RG equation (16) still holds. As a consequence, the critical exponents below d u are the 
same than for m = 2. The steady-state concentration can be written (for d < d u ) as: 

/ Jx d/(d+2) 

c to ,r(oo) = A m l — J , 

and the universal amplitude A m can be computed within an e-expansion. The approach 
to the steady-state is still exponential with fc-dependent amplitudes. Extension to the 
reactions mA — > I A with / < m (m > 2) is also possible (see 0]). 

Another natural extension of this work is to consider the possibility of reversible 
reaction (for example A + A # A). For such system a completely different physics is 
expected. Indeed, a quick investigation of the associated action shows that the upper 
critical dimension is no longer 2 but 4. In addition, it appears that a wavefunction 
renormalization is needed, giving rise to anomalous dimensions. The computation of 
scaling functions for this problem is currently under investigation. 



Appendix A. 

We aim to show that in the limit J — > oo, the e-expansion breaks down. For this purpose, 
let us consider the action (|T8| ) obtained by the shift of the field a by its classical value. 
The shift permits us to suppress the source term in the action; in other words, shifting 
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the field is equivalent to perform the infinite sum generated by the source term. One 
could then think that any quantity can be calculated for arbitrary J. This is however 
not true for the following reason. To calculate a given quantity, we expand the action 
with respect to the \ff{a 2 cl + 2a c \7] + rj 2 ) term. We shall then obtain a power expansion 
in term of A with coefficient proportional to a 2 cl . But for large J, a c \ tends to (J/A) 1 / 2 . 
The expansion in A is partly replaced by an expansion in J, which does not lead to an 
e-expansion. The only possibility to revert this would be to treat non-perturbatively 
the action, which is out of question. 



Appendix B. 

The problem is to compute the e-expansion of 

1 dxx~ 1+e f(x). (Bl) 



o 



Putting naively x~ 1+e = + e lnx + 0(e 2 )] in flET| ) clearly fails, because if /(0) ^ 0, 
Jq dx x~ 1 f(x) diverges. One way to avoid this problem, is to treat x~~ 1+e as a generalized 
function. We shall not give here a detailed discussion of the generalized functions 
(see |20| for an introduction). We only quote the result, namely: 

f'dxx-^fix) = -/(0) + f'dxx^ifix) - /(0)] 
Jo e jo 

+e f dxx'Hnxlfix) - f(0)} +0(e 2 ). 
Jo 

Remark that there is a pole of order 1/e. 
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